Forum for Electromagnetic Research Methods and Application Technologies (FERMAT) 


Review of Modified Algorithms for Synthetic 
Aperture Radar Imaging at High Squint Angles 


Po-Chih Chen and Jean-Fu Kiang 


Abstract—This paper summarizes the modified versions of 
chirp scaling algorithm (CSA), range-Doppler algorithm (RDA) 
and frequency-domain algorithm (FDA) for synthetic aperture 
radar (SAR) imaging at high squint angles. The reviewed 
methods can reconstruct better images with less computational 
load and memory than their conventional counterparts, which is 
verified by simulations. SAR imaging with highly-elliptical-orbit 
(HEO) satellites based on these techniques is also reviewed. 


Index Terms—synthetic aperture radar (SAR), high squint 
angle, chirp scaling algorithm (CSA), range-Doppler algorithm 
(RDA), frequency-domain algorithm (FDA), highly-elliptical- 
orbit (HEO). 


I. INTRODUCTION 


Synthetic aperture radar (SAR) techniques have been widely 
used to acquire high-resolution images of the ground. Many 
conventional SAR imaging algorithms work properly when the 
squint angle is small. However, as the squint angle becomes 
large, the slant range becomes a more complicated function of 
relevant geometrical parameters and the range-azimuth cou- 
pling also grows. To deal with such issues, more complicated 
algorithms were proposed, such as azimuth nonlinear chirp 
scaling (ANCS) [1], [2], frequency nonlinear chirp scaling 
(FNCS) [3], fractional CSA (FrCSA) [4], modified CSA [5], 
[6], range migration algorithm (RMA) [7], [8], improved step 
transform (IST) [9], azimuth overlapped subaperture algorithm 
(AOSA) [10], and so on. 

At high squint angles, the range-Doppler spectrum resides 
on a very oblique subdomain and may fold in the Doppler- 
frequency dimension. A higher pulse repetition frequency 
(PRF) is required to remedy this issue [11], [12], which 
increases computational load and memory size. In addition, 
the difference between the maximum and the minimum slant 
ranges to the target area becomes very large, which demands 
larger matrices to store the received signals and the interme- 
diate data, respectively [9]. In [13] and [14], modified SAR 
imaging algorithms were proposed to reduce the memory size 
by rotating the received baseband signals in the time domain to 
a new coordinate system in which the slant ranges at different 
slow times fall in similar ranges. The corresponding spectrum 
is also rotated to reside on a nearly rectangular subdomain, 
and the folded-spectrum problem is also resolved. 

Inspired by these modified SAR imaging algorithms, highly- 
elliptical-orbit (HEO) satellites [15], [16], conventionally used 
for communications over high-latitude areas, are envisioned 
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as platforms for SAR imaging. In [17], a constellation of two 
satellites in the Tundra orbit were simulated for SAR imaging. 
The Doppler centroid is extremely large when the satellite 
moves far away from the apogee, due to high eccentricity of 
the orbit, demanding a large number of samples along the 
slow-time dimension. In addition, the difference between the 
maximum and the minimum distances between the satellite 
and the target area is extremely large, demanding a much larger 
number of samples in the range dimension. Similar to the high- 
squint SAR scenarios, a modified frequency-domain algorithm 
(FDA) was proposed for the HEO SAR missions, which takes 
much smaller computational load and memory size than the 
conventional FDA. 

This review paper is organized as follows. Four CSA-based 
methods [13], designed for spaceborne SAR imaging at high 
squint angles, are reviewed in Sections II; an RDA-based 
method [14] is reviewed in Sections III; an FDA-based method 
[17], tailored for HEO SAR imaging, is reviewed in Sections 
IV. Finally, a brief summary is presented in Section V. 


II. CSA-BASED METHODS 








(0, Yb, h) 


(Xe, Ye, 0) 


Fig. 1. Flight path of a platform carrying a SAR radar that points at a look 
angle 0p and a squint angle 0s towards the target area [13]. 


Fig.l shows the flight path of a platform carrying a SAR 
radar that points at a look angle 6, and a squint angle 6, 
towards the target area [13]. The platform flies in the y 
direction from (0, y», h) to (0, ye, h), at velocity of Vp and 
height of h. The beam center point (BCP) falls at (£e, Ye, 0), 
with ze = htanĝe and ye = htand,/cos6¢. The slant 
range from the center of the flight path to the BCP is Rso = 
\/x2 + y2 + h?, and the range to the BCP is Ro = \/x?2 + h?. 
At azimuth time 77, the slant range from the platform to a point 
target at (x, y,0) is Re(n,2,y) = fz? + (nVp — y)? +h?, 
and the range to the point target is R(x) = Va? + h?. 
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Fig. 2. Flowcharts of CSA-based methods for high-squint SAR imaging, (a) 
method A and (b) method B [13]. The processes enclosed with parentheses 
are different from those in the conventional CSA. 


Fig.2(a) shows the flowchart of CSA, method A for high- 
squint SAR imaging [13], and the processes different from 
those in the conventional CSA are enclosed with parentheses 
[18]. The compensation of coupling is essential to high-squint 
SAR imaging, which is implemented in the two-dimensional 
frequency domain, requiring additional range FFT and range 
IFFT. 

The transmitted linear frequency-modulation (LFM) signal 
is represented as s,(T) = we(r)el2"foTt+I™ rr” where fo 
is the carrier frequency, AK, is the FM rate, 7 is the range 
time, we(T) = rect(7/T;) is the range envelope with duration 
T;, and rect(7) is a rectangular function, which has unity 
amplitude in |r| < 1/2 and zero otherwise. The demodulated 
scattering signal, at range time 7 and azimuth time 7, from a 
point target at (x, y,0) is expressed as 


Srb ce n) = Ao(a, y)We(T = 2R Ag, x, y)/C) 
—jAnfoRs(n,x,y)/cetjnK,[rT—-2Rs(n,x,y)/c]? 


(1) 


where Ao(x, y) is a complex amplitude, which is arbitrarily 


€ 


set to one for convenience, and the radiation pattern from the 
SAR radar towards the point target is approximated as one. 

Next, take the Fourier transform of the baseband signal in 
(1) with respect to 7 and 7) to obtain 


Se (lees fr) = C1C2 Ag (x, y)Wel fre ati / Kr e520 Fn (u/ Ve) 





AnR(a oie 3 
C fo Jó 
where W.(f-) = we(f-/Kr), ci and cg are constants 
Cis . j 
of integration, and D = 1 — IVA fE Typically, D4 > 
of, fP 
Fo + 7 , thus the square-root term in (2) can be approxi- 
0 
mated as : 
Ol i Oe. Ie 
D445 2¢D4+ 0 ~  t F®) 
fo off Dfo 8D°V2 fo  8D°VE fe 


As a comparison, only up to f? terms are used in the 
conventional CSA. By substituting (3) into (2), we have 


(4) 


To compensate for the f?-term in the exponent of (4), a 
coupling compensation filter is devised as 


tok fe 
Helden) = exp poit 


So C= tn) a C1C2 Áo (2; y)We(frje attr / Kr e~ 42 Fn (u/ Vo) 


aod GEOD |p, f Ehh | int: 
ri Dfo  8D8V2fe ° 8D>V2 fe 








which can be ignored in low-squint SAR imaging. The signal 
after applying the coupling compensation filter is inverse 
Fourier transformed with respect to f,- to obtain 


Sa(r, Fn) = €1€2C3A0(2, yje 7I?" nly Ve) 


u (ez h- ZR] ) ow {S| 





2R(x)]? 
exp jak, z- =| i (5) 
cD 
here cs i tant of integration, K,, = os d 
W cz 1S a constant of integration, K,, = I KZ an 
_ Ro) fp 
2D3V?2 fe 


The slant range in the range-Doppler domain is decomposed 
as r.(R(x), fy) = R(x)/D(f,), and the total RCM of a 
point target is defined as RCM;(R(x), fy) = rs(R(x), fy) — 
rs(R(x), fy,.), where f,,. is commonly chosen at the Doppler 
centroid. The total RCM is further decomposed into a bulk 
RCM, RCM, (f;,) = RCM: (Rr, fn), determined at a reference 
range R,, which is commonly chosen to be the range to 
the BCP, Ro = \/x2+h?; as well as a differential RCM, 
RCMag(R(2), fy) = RCM:(R(x), fy) — RCMo(f,). The 
range time to the point target, after compensating the differen- 


tial RCM, becomes T = F [rs (R(x), fn) — RCMa(R(zx), fn). 
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2 
The range time to the reference point is T, = — [rs(Ry, fy). 
The time shift associated with the differential RCM is Ar = 


- x RCMa(R(2), fn) = Po 





—1|7’, where 7’ = 7 — 


T — 
cD( fn) 
To compensate for the differential RCM, a chirp scaling 
function is devised as 





= 


Hales Jn) = exp laf fsc(u, fau} 


where fsc(u, fn) = KmAT(u, fi). If Km is approximated as 
independent of range, the scaling function can be reduced to 


Hc (T, fn) = exp {in ie E r ( p a 


which is multiplied with (5) to derive a scaled signal in the 
range-Doppler domain as 





S5(T, fn) = C1 C2C3 Áo (T, yje 7I? fn U/W) 
1 2R(x) —j4r R(x) foD/c 
We (7 z DE |) e 


a {ist -20 | 


exp {in ee = 1 (- — =) (6) 


Next, by taking the Fourier transform of the scaled signal 
in (6) with respect to 7, and imposing the approximations that 
Z~0 and fy, ~ fn, we obtain 
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which is multiplied with a range-compensation filter 
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to derive 
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which is then inverse Fourier transformed with respect to f- 
to have 


Sg (T Ja = C1 C2C3 Áo (x, ye J27 Fn (y/Vp) p—J47 R(x) foD/c 
e AnKm |R) Rr ° 1 D 
x i = 
E P 


F,sinc |r, z 7 are 


with —F,/2 < fr < F,./2 — F,./N,, where F;, is the range 
sampling rate. 
By applying an azimuth-compensation filter 


Hac( fn) = exp joer 











(8) 


and a residual-phase compensation filter 








Hype( fin) = eXD {ee = - Z| | i - i. 


(9) 
to Ss(T, fn), we obtain 
So(r, Fn) = €1€2C3A0 (2, ye F27 Fn ule) 
cana 


which is inverse Fourier transformed with respect to fy to 
obtain 


s10(T, 1) = c1c2c3Ao (z, y)F,sinc fr, | E AL ) 


cD( fo) 
Fasinc {Fa (Nn — y/Vp)} (10) 


with —Fa/2 < fn < Fa/2— Fa/Na, where Fa is the azimuth 
sampling rate. Based on (10), the 7-7) coordinates of the point 


2R 
target are estimated as T = 2R and n = E which are 
cD( fin, ) Vp 
ED( fr )T i 
mapped to £ = eT —h? and y = Vpn 


B. Method B 


Fig.2(b) shows the flowchart of CSA, method B for high- 
squint SAR imaging. The processes before chirp scaling 
remain the same as in method A. A two-dimensional similarity 
theorem is applied to reduce computational load and memory 
size, which states that if the signals in the 7-7 plane is ro- 
tated by a given angle, the corresponding Fourier-transformed 
signals in the f--f, plane will be rotated by the same angle. 

The received baseband signals, represented in (1), are stored 
in a matrix of dimension NV, x N,, where the azimuth sam- 
pling number (Na) determines the resolution in the azimuthal 
direction, and the range sampling number (N,.) must be large 
enough to store all the scattered signals from the target area. 
For the convenience of applying the FFT algorithm, both N, 
and N, are incremented, respectively, to the nearest integers 
that are powers of two. 

Figs.3(a) and 3(b) show the storage areas of received 
baseband signals in (7 — To, n) plane and (r’ — To, n’) plane, 
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Fig. 3. Storage area of received baseband signals in (a) (7 — To, 7) plane 
and (b) (7’ — To, n’) plane. 


respectively, where 7 = 2R,0/c. The number of range 
samples in Fig.3(b) is much smaller than that in Fig.3(a). 
The rotation angle 6, between these two coordinate systems 
4 20/6 

L/Vp 
length of the flight path, and D = \/x2 4+ (ye — yo)? +h? — 

£2 + (Ye — Ye)? + h? is the length difference between the 
two line segments that connect the BCP to the starting point 
and ending point, respectively, of the flight path. Explicitly, 
the (T — 7,7) coordinates are related to the (T° — 70,7’) 
coordinates as 


is chosen as 6, = tan` 





, where L = ye — yp 1S the total 


T (T — T)) cosO,. — 7’ sin 6, + To 
n (T — To) sin 6, +7’ cos 6, 
g(t’, 17’) 
h(a’, n’) | 


By taking an inverse Fourier transform of the scaled signal 
in (6) with respect to f,, we have 


ST. n) = T {55 (7, In)} 


The scaled signals in the (7,7) coordinates can now be 
represented in the (7’, 7’) coordinates as 


(11) 


Ssc (es n’) = Ssc (T, n) 
T=9(7'.n'),n=h(1' n’) 


TABLE I 
PARAMETERS OF SAR MISSION [18]. 
effective radar velocity V; 
chirp pulse duration 
bandwidth 
range sampling rate 


mk 


: 


* a/b: a is used at 0s = 60° and 6 is used at 0; = 80°. 


10.42 / 41.67* | ns 


6,800 7 1, 700° 
0.1570.59* | ms 
16,384 





FERT 


which are stored in a matrix of dimension N} x N/.. Without 
loss of generality, we choose Ar’ = Ar and An’ = An. 

The filters Hac(fn) in (8) and H,pc(fn) in (9) can be 
relabeled as Hac( fr, fn) and Hrpc(f-, fy), respectively. Both 
filters and the range-compensation filter in (7) can be repre- 
sented in the f,-f,, plane as 


Weal ies fo) = Hrc( fr, fn) 





fr=91 (ft fh) fn=hi Falan 





Hal -s fa) = Hach fr, fn) 
fr=9il ft fn) fn=hi (ft Sa) 


Hse es a = ype fr, tn) 





fr=9i (ft) fn=hi (ft fp) 


where (f+, fn — fac) are related to ( f}, fp — fac) by the same 
rotational matrix in (11) as 


f {cos 0, — (f%, — fac) sin, 


Ín f, sin 6, + (ff, — fac) cos Or + fac 
Cita) 


hı (fy, fa) 


C. Methods C and D 


Method C is revised from method A by implementing 
the range-compensation process right after the coupling- 
compensation process, eliminating a range FFT and a range 
IFFT processes. Method D inherits the advantage of reduced 
memory storage in method B and implements the range- 
compensation process right after the coupling-compensation 
process, as in method C. In addition, the coupling compen- 
sation filter is implemented in the rotated coordinates, further 
reducing the computational load. 

Two missions with squint angles of 60° and 80°, respec- 
tively, are simulated to verify the efficacy of these methods. 


D. Mission with 0, = 60° 


Table I lists the parameters of SAR missions simulated in 
this work [18]. Fig.4(a) shows the received signals, computed 
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Fig. 4. (a) Received signals in 7-7 plane and (b) image reconstructed with 
method A, 0s = 60°, Nr = 16,384, Na = 16,384 [13]. 
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TABLE II 
PERFORMANCE INDICES OF SAR IMAGING, 0, = 60° [13]. 





by using (1), in the 7-7 plane. Fig.4(b) shows the image 
reconstructed with method A, in which local maxima are 
observed along the x axis and a y’ axis. 

Fig.5(a) shows the scaled signals in the r’-ņ' plane, and 
Fig.5(b) shows the image reconstructed with method B. The 
number of range samples is reduced to N/ = 4,096, and the 
number of azimuth samples remains as NV’ = 16,384. The 
image reconstructed with method B appears similar to that 
with method A. The memory required by method A to store the 
received signals is on the order of N, x Na, while that required 
by method B is on the order of N} x N}. The reconstructed 
images with methods C and D also appear similar to those 
obtained with methods A and B. 

Table II lists the performance indices of SAR imaging by 
using different methods, with 6, = 60°. The impulse response 
width (IRW) is defined as the separation between two half- 
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Fig. 5. (a) Scaled signals in 7’-7’ plane and (b) image reconstructed with 
method B, 6, = 60°, N/ = 4,096, N/, = 16,384 [13]. 
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power points, one on each side of the peak-intensity point; 
the peak sidelobe ratio (PSLR) is the ratio between the peak 
intensity and that of the strongest sidelobe; the integrated 
sidelobe ratio (ISLR) is the ratio between the total intensity 
in the main-lobe and that in all the sidelobes; and the target 
offset (TO) is the offset of the peak-intensity point in the 
reconstructed image from the actual point target. The images 
reconstructed with these four methods have similar IRW, 
PSLR, ISLR and TO, in both x and y’ directions. However, 
methods B and D take only 25 % of memory as compared 
to methods A and C, due to the signal-rotation scheme. The 
numbers of multiplications (NOMs) required by methods A, 
B, C and D have the ratio of 1 : 0.884 : 0.687 : 0.568. In 
summary, method D out-performs the other three in NOMs 
and memory. 


E. Mission with 0, = 80° 


Fig.6 shows the image reconstructed by using method A, 
with 86, = 80°. When method B is applied, the number 
of range samples is reduced to Ni = 1,024, and that of 
azimuth samples remains as N’ = 16,384. Table III lists 
the performance indices of SAR imaging by using different 
methods, with 0, = 80°. The IRWs in the x direction (IRW,,) 
with methods A and B are smaller than those with methods C 
and D. The images reconstructed with these four methods have 
similar IRW,,,, PSLR,, ISLR,, TO; and TO. However, 
methods B and D take only 6.25 % of memory as compared 
to methods A and C, due to the signal-rotation scheme. The 
NOMs required by methods A, B, C and D have the ratio 
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Fig. 6. Image reconstructed with method A, 0, = 80°, Nr = 16,384, 
Na = 16, 384 [13]. 


TABLE III 
PERFORMANCE INDICES OF SAR IMAGING, 6, = 80° [13]. 


1.5876 1.0568 1.5876 





of 1 : 0.708 : 0.687 : 0.387. In summary, method D out- 
performs the other three in NOMs and memory, but is poorer 
than methods A and B in IRW,,. 


III. RDA-BASED METHOD 
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Fig. 7. Flowchart of RDA-based method for high-squint SAR imaging 
[14]. The processes enclosed with parentheses are different from those in 
the conventional RDA [18]. 


Fig.7 shows the flowchart of an RDA-based method for 
high-squint SAR imaging [14]. To begin with, the time-domain 
received baseband signals are rotated to a new coordinate 
system. Two major processes, range compression and azimuth 
compression, are implemented in the f,-f, domain. Four 
relevant filters are applied in the two-dimensional frequency 
domain after rotation. At last, the processed signals are rotated 
back to the 7-7) domain and mapped to the x-y coordinates. 


The received baseband signals stored in the (7,7) coordi- 
nates are transformed to the (7’, 7’) coordinates as 


(12) 
t=g(7' sn"), n=h(7',n') 


Srp(7', 7) = Srb(T, 7) 


which are interpolated to a uniform grid in the r’-7’ plane 
and stored in a matrix of dimension NV} x N. Without loss 
of generality, we choose Ar’ = Ar and An’ = An. The 
range-compression filter used in the conventional RDA is [19] 


Aa) = lth /Km 


where 


1 AoR(x) fF 


K Dy 


l — 
Km 





which is applied in the fr-fn domain instead of the fr-n 
domain as in the conventional RDA. An azimuth-compression 
filter is devised as [19] 


Hacl fr, fn) = clam R(a)D/ro 


which is applied in the /,-f, domain instead of the 7-f,, 
domain as in the conventional RDA. A range cell migration 
correction (RCMC) filter is devised as [19] 


Hreme(fr, fn) = exp l; mon > p ay 


At high squint angles, a coupling-compensation filter is re- 
quired [19] 


C 


TrR(x) f? ff 


HAT, Ja) = exp 2D5 fS V2 
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These four filters are rotated to the f,-f,, plane as 


Fits a = Hye (fr; fn) 





f-=91 (ft fh) fn=hi aa) 
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Toome J ~ 7) — Pisce: To) 





fr=91 (ft fn) fn=hi (testa) 
Hal fa) = Hach fr, fn) 





fr=91 (F157), fn=h1 FSi) 
Two missions with squint angles of 60° and 80°, respec- 
tively, are simulated to verify the efficacy of this method. 


A. Mission with 0, = 60° 


The same parameters of SAR missions, as listed in Table I, 
are used in the simulations [18]. Fig.8(a) shows the received 
signal in the 7-77 plane, computed by using (1). Fig.8(b) shows 
the image reconstructed with the conventional RDA, in which 
local maxima are observed along the x axis and a y’ axis. 

Fig.9(a) shows the received signal in the 7’-7)’ plane, com- 
puted by using (12); and Fig.9(b) shows the image recon- 
structed with the RDA-based method. The number of range 
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(a) Received signal in T-7) plane and (b) image reconstructed with 
Nr = 16,384, Ng = 16, 384 [14]. 


Fig. 8. 
the conventional RDA, 6, = 60°, 


TABLE IV 
PERFORMANCE INDICES OF IMAGING SAR, 0, = 60° [14]. 


IRW; (m) 9.6928 9.6833 9.6826 
PSLRz (dB) —13.2521 —13.2418 —13.2406 
IRW, (m) 16.3341 16.1536 8.1589 


PSLR,,, (dB) —13.3216 —13.2337 —13.1988 
NOM (million) | 35,165 11,542 23,622 
memory (GB) 





samples is reduced to N/ = 4, 096, and the number of azimuth 
samples remains as N’ = 16,384. This scenario is labeled 
as scheme A. The image reconstructed with the RDA-based 
method appears similar to that with the conventional RDA. 
The memory required by the conventional RDA and the RDA- 
based method to store the received signals are on the order of 
N, x Na and N} x N}, respectively. 

A second set of parameters, NV’ = 4,096 and N/ = 32, 768, 
are chosen to form scheme B. The number of range samples is 
the same as that in scheme A, but that of the azimuth samples 
is doubled. The NOMs and memory required in scheme B 
are larger than that in scheme A. Fig.10 shows the image 
reconstructed with the RDA-based method. The number of 
range samples required is only a quarter that required by the 
conventional RDA, and the number of azimuth samples is 
twice that of the conventional RDA. The reconstructed image 
appears to have finer spatial resolution in the y’ direction 
because the synthetic aperture length is doubled. 

Table IV lists the performance indices of SAR imaging by 
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Fig. 9. (a) Received signal in 7’-ņ’ plane and (b) image reconstructed with 


the RDA-based method in scheme A, 6s = 60°, N} = 4,096, N} = 16,384 
[14]. 
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Fig. 10. Image reconstructed with RDA-based method, scheme B, 0s = 60°, 
Ny = 4,096, NZ = 32, 768 [14]. 


using the RDA-based method, with 6, = 60°. The images 
reconstructed with the conventional RDA and the RDA-based 
method in scheme A are similar. However, the latter takes 
only 33 % of NOMs and 25 % of memory as compared to 
the former. The spatial resolutions in the x direction with the 
conventional RDA and the RDA-based method are similar in 
both schemes, and the spatial resolution in the y’ direction with 
the latter in scheme B is only half that with the former. The 
RDA-based method in scheme B takes 205 % of NOMs and 
twice memory as compared to those in scheme A. However, 
the former takes only 67 % of NOMs and one half memory 
as compared to the conventional RDA. In summary, the RDA- 
based method in scheme B out-performs the conventional RDA 
in spatial resolution along 4’ direction, NOMs and memory. 
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B. Mission with 6, = 80° 
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Fig. 11. Image reconstruction (a) with conventional RDA, 0s = 80°, Ny = 
16,384, Na = 16,384; (b) with RDA-based method in scheme C, 0s = 80°, 
Ni = 1,024, N! = 16,384; (c) with RDA-based method in scheme D, 
0s = 80°, N}. = 1,024, N} = 65,536 [14]. 


Next, consider a scenario with the squit angle increased to 
0s = 80°. Fig.11(a) shows the image reconstructed with the 
conventional RDA. Fig.11(b) shows the image reconstructed 
with the RDA-based method in scheme C, with N’ = 1,024 
and N} = 16,384. The number of azimuth samples in scheme 
C remains the same as N} = Na, and that of range samples 
is only 1/16 that of the conventional RDA. 

In scheme D, N} = 1,024 and N? = 65,536, the number 
of range samples remains the same as that in scheme C, 
while that of azimuth samples is increased by four folds. 
Fig.11(c) shows that the reconstructed image appears to have 
finer spatial resolution in the y’ direction because the synthetic 
aperture length is increased by four times. 

Table V lists the performance indices of SAR imaging in 


TABLE V 
PERFORMANCE INDICES OF SAR IMAGING, 06, = 80° [14]. 


IRW; (m) 17:3193 17.3099 17.4109 
IRW, m) 66.3484 66.4359 16.3679 


PSLR,,, (dB) —13.2816 —13.2475 —13.1884 
NOM (million) | 35,165 2,751 11,542 
memory (GB) 





schemes C and D, with 6, = 80°. The image reconstructed in 
scheme C appears similar to that with the conventional RDA. 
However, the former takes only 7.8 % of NOMs and 1/16 
of memory as compared to the latter. The spatial resolutions 
in the x direction with the conventional RDA and scheme D 
are similar, and the spatial resolution in the y’ direction with 
the latter is only half that with the former. Scheme D takes 
419 % of NOMs and four times of memory as compared to 
scheme C. However, the former takes only 33 % of NOMs 
and one quarter of memory as compared to the conventional 
RDA. In summary, the RDA-based method in scheme D out- 
performs the conventional RDA in the spatial resolution along 
y’ direction, NOMs and memory. 


IV. FDA-BASED METHOD 


received baseband signal 
(signal rotation) 
2D FFT 


( modified ) 


RFM filter 


final image 





Fig. 12. Flowchart of FDA-based method for high-squint SAR imaging 
[17]. The processes enclosed with parentheses are different from those in the 
conventional FDA [20]. 


Fig.12 shows the flowchart of the FDA-based method for 
high-squint SAR imaging [17], and the processes different 
from those in the conventional FDA [20] are enclosed with 
parentheses. The received baseband signal is rotated in the 7- 
7 plane before being Fourier transformed to the fr-fn plane. 
The reference function multiply (RFM) filter is rotated in the 
f--fn plane by the same angle. Finally, the processed signal 
is reversely rotated in the 7-7 plane to reconstruct the final 
image. 

The slant range in a SAR mission can be represented as 
[20] 


(13) 


Forum for Electromagnetic Research Methods and Application Technologies (FERMAT) 


where 7;(7) and 7,(7) are the positions of the satellite and 
the target, respectively, at azimuth time n. Define 


R(n) = Fs (n) — Fa (n) 

V(n) = a = ts() — Bg(n) 

A(n) = a = Gs() — Gg(n) 

B(n) = A = bs (n) — b(n) 

D(n) = ana = ds(n) — dg(7) 
dD(n) 





where U,(7) and (n) are the velocity vectors of the satellite 
and the target, respectively; @,(7) and @ (7) are the accel- 
eration vectors of the satellite and the target, respectively; 
b,(7) and b,(7) are the time derivatives of a;(7) and G,(n), 


respectively; ds(7) and d,(n) are the time derivatives of 


b,(7) and bg(7), respectively; €,(7) and €,(7) are the time 
derivatives of d,(7) and d,(m), respectively; all at azimuth 
time n. The slant range can be represented in terms of a fifth- 


order Doppler parameter range model (DRM-5) as [20] 


= Re+ X knn” (14) 
where R, = | R(0)| is the slant range at 7 = 0, and 
_ R(0)-V(0) 
"= RO oe 
Le R* (0) - A(0) + V* (0) - V(0) — ky 
i 2|R(0)| 
ia R'(0) - B(0) + 3V*(0) - A(0) — 6ky ke 
6|R(0) 
Te R+ (0) - D(0) + 4V* (0) - B(0) + 34+ (0) - A(0) 
a 24| R(0) 
k k- + 2kık3 
2|R(0)| 
Bec R'(0) - E(0) +5V*(0) - D(0) + 10.4*(0) - B(0) 
T 120| R(0)| 
_ kok3 + kika 
|R(0)| 


An RFM filter is devised as [20] 


Are (f-, Ty) = oe Jlér(fr)+¢Rrem (fr >fn)+ba(fn)+be( fr fndI(1 5) 


where $,(f7), Orom (fr, fn)» Galfn) and (fr, fn) are the 


phases to implement range compression, range cell migration 
(RCM) correction, azimuth compression and decoupling, re- 


TABLE VI 
PARAMETERS OF TUNDRA ORBIT [16]. 


semi-major axs — |a _| 42,164 | km | 
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The RFM filter is applied in the f}-f, domain, with the form 
of 


Hrrm( fr fa) (16) 


fr=91 (F4, fh), fn=hi af 


= Hrrm( fr, fn) 


A. Satellite Orbit and Ephemeris 


Table VI lists the parameters of a Tundra orbit [16]. A 
constellation of two satellites are deployed to detect a point 
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target specified in the Earth-centered Earth-fixed (ECEF) co- 
ordinates. The right-of-ascending-node and the true anomaly 
of the second satellite differ from those of the first one by 
180°, and the other parameters are the same. 
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Fig. 13. (a) Ephemeris of two Tundra satellites, the solid circle marks a point 
target. (b) Time history of squint angles measured from two Tundra satellites 
[17]. 


Fig.13(a) shows the ephemeris of the two Tundra satellites. 
The solid circle marks a point target selected for simulation. 
Along each ephemeris, the thick part marks the working 
window, which is centered around the apogee and lasts for 
about half of the period. Fig.13(b) shows the time history of 
squint angles (@,) measured from each of the two Tundra 
satellites. The period of each Tundra satellite is 7’ = 24 
hours, and the first satellite passes the apogee at t = 0. In 
other words, the first satellite operates in 0 < t < T/4 and 
3T/4 <t < J, and the second one operates in the rest of the 
period. 


The SAR images of the target point are taken at six instants, 
marked as events /\-&. in Fig.13. The four events F}- 
E, involve high squint angles, event E5 involves a small 
squint angle, and event fg involves zero squint angle. The 
conventional FDA is also applied to reconstruct images for 
comparison. Events Ł1-E4 impose more challenges than the 
other two events on applying the HEO satellites for SAR 
imaging. 


TABLE VII 
PARAMETERS OF GEO SAR MISSIONS [20]. 
chirp pulse duration 
bandwidth 


a/b*: a is used in Events E;-Ey4 and b is used in Events Es, Ee. 





B. Event Ey 


Table VII lists the parameters of a GEO SAR mission [20], 
which are used for the HEO SAR simulations. The memories 
required by the conventional FDA and the FDA-based method 
to store the received signal are on the order of N, x Na and 
N’. x N}, respectively. 
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Fig. 14. (a) Received signal in 7-7) plane and (b) image reconstructed with 
conventional FDA, Nr = 16,384, Na = 8,192 [17]. 


Fig.14(a) shows the received signal in the 7-7 plane, which 
is computed by using (1). Fig.14(b) shows the image recon- 
Structed with the conventional FDA, in which local maxima 
are observed along the x axis and a y’ axis. 

Fig.15(a) shows the image reconstructed with the FDA- 
based method. Scheme A is formed by reducing the number 
of range samples from N,. = 16,384 to N, = 4,096, while 
keeping the number of azimuth samples as N’ = 8,192. 
The reconstructed image appears similar to that with the 
conventional FDA. 
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Fig. 15. Image reconstructed with FDA-based method, (a) scheme A, Nj. 
4,096, Ni, = 8,192, (b) scheme B, N/. = 4,096, N/ = 16,384 [17]. 


TABLE VIII 
PERFORMANCE INDICES OF SAR IMAGINGIN EVENT £1 [17]. 


index conventional | scheme A scheme B 


IRW; (m) 29.651 29.9558 30.6552 
IRW, (m) 220.0219 222.4984 112.075 


PSLR,,, (dB) —13.2704 —13.2437 —13.2968 
NOM (million) 15,166 4,932 10,133 
memory (GB) 





Scheme B is formed by choosing N/ = 4,096 and in- 
creasing N’ from 8,192 to 16,384. Thus, the NOMs and 
memory required in scheme B are larger than those in scheme 
A. Fig.15(b) shows the image reconstructed with scheme B. 
The number of range samples is only a quarter that of the 
conventional FDA, and the number of azimuth samples is 
twice that of the latter. The reconstructed image has finer 
spatial resolution in the y’ direction because the synthetic 
aperture length is doubled. 

Table VII lists the performance indices of SAR imaging 
in event /,. The images reconstructed with the conventional 
FDA and the FDA-based method in scheme A appear similar. 
However, the latter takes only 33 % of NOMs and 25 % of 
memory as compared to the former. The spatial resolutions in 
the x direction with both methods are similar; and the spatial 
resolution in the y’ direction with the FDA-based method 
in scheme B is about half that with the conventional FDA. 
Scheme B takes 205 % of NOMs and twice memory as 
compared to scheme A. However, the former takes only 67 
% of NOMs and one half of memory as compared to the con- 
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TABLE IX 
PERFORMANCE INDICES OF SAR IMAGING IN EVENT E; [17]. 


FDA based 
IRW: m) 33.425 33.2845 


PSLRz (dB) —13.0985 —13.1341 
IRW,,, (m) 144.074 145.2704 


PSLR,,, (dB) —13.2202 —13.2274 
NOM (million) | 31,407 20,803 
memory GB) [4 P 





ventional FDA. In summary, the FDA-based method in scheme 
B out-performs the conventional FDA in the spatial resolution 
along y’ direction, NOMs and memory. The simulations in 
events H5-H, demonstrate similar results. 


C. Event Es 
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Fig. 16. Image reconstructed (a) with conventional FDA, Ny = 16,384, 
Na = 16,384, (b) with FDA-based method, Ne = &, 192, N? = 16,384 


[17]. 


Fig.16(a) shows the image reconstructed with the conven- 
tional FDA, in which local maxima are observed along x and 
y’ axes. Fig.16(b) shows the image reconstructed with the 
FDA-based method. The number of range samples is reduced 
from N, = 16,384 to 8,192, while the number of azimuth 
samples remains as N} = 16,384. The image reconstructed 
with both methods have similar quality. 

Table IX lists the performance indices of SAR imaging in 
event £5. The FDA-based method takes only 66 % of NOMs 
and one half memory as compared to the conventional FDA. 
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V. SUMMARY 


We have reviewed four modified chirp scaling algorithms 
(CSAs), one modified range-Doppler algorithm (RDA), and 
one modified frequency-domain algorithm (FDA), which are 
used to improve the quality of SAR imaging under high squint 
angles. Scenarios with squint angles up to 60° and 80° have 
been simulated to demonstrate the efficacy of these methods 
in terms of spatial resolution, computational load and memory 
size. Two HEO satellites in a Tundra orbit are envisioned 
for SAR imaging, which are capable of providing a complete 
coverage over high-latitude areas. 
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